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The properties of a strongly interacting quark plasma are investigated by molecular 
q ; dynamics method including non-abelian quark-quark potential. Our main goal is to study 
the thermalization process in this system. We find an interesting resonance-like behaviour: 
at a characteristic time close to the inverse plasma frequency the quark plasma is heated 
up substantially via energy transfer from quark potential energy into one particle kinetic 
energy. Color rotation mechanism enhances the effectivity of this heating process, leading 
^ ' to a very fast thermalization with high temperature. 



Introduction 



Early theoretical speculations suggested the appearance of a weakly interacting per- 
turbative quark-gluon plasma state in heavy ion collisions at ultrarelativistic energies. 
However, recent data obtained at RHIC experiments, and their theoretical analysis, indi- 
cate the formation of a strongly interacting deconfined matter [[!]. In fact, lattice-QCD 
calculations on equation of state have already shown the strong coupling and the the- 
oretical analysis based on massive quasi-particle picture described quantitatively many 
properties of this matter Thus we expect a relatively dilute, massive, quark and 

antiquark dominated matter to appear before hadronization. Widely used molecular dy- 
namics simulation [El El is the ideal tool to study such a particle compound, to investigate 
non-perturbative features and to determine different properties in the strongly interacting 
matter. In this simulation we can study the influence of the color charge fluctuation and 
pair correlation function can be generated, displaying if the quark matter is gaseous or 
rather liquid like. Although molecular dynamics models for non-abelian interactions are 
very much approximate, we can study qualitatively very important microscopic processes, 
for what other methods do not exist. Here we summarize the main findings of our work, 
details can be found in Ref. 
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2. The Molecular Dynamics Model 



In the model the quarks are the only actively simulated species, the gluons are rep- 
resented as a background field. This approach is similar to the case of the classical 
multicomponent or one component plasma (OCP) model, where the motion of the ions is 
of interest, the electronic background is treated as a constant field [HI Ell- 
in the model pairwise interaction of the quarks is assumed. The possible realizations 
of (single flavored) two-quark (QQ) system are composed of symmetric sextet and anti- 
symmetric anti-triplet combinations. 

The quark-quark effective interaction is repre- 
sented by a Coulomb-like potential in the form of 



Table 1: (Aj • Xj) values for QQ-pairs. 
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acting between the i th and j th particle, where a 
g 2 /47r is the interaction coupling constant, is tl 
inter-particle distance and (A, • Xj) depends on the color states of the particles as shown 
in Table [U The parameter D introduced in Table [T] represents values of (A, ■ Xj) = +| or 
— | randomly chosen with equal probabilities, due to the random realization of symmetric 
and antisymmetric QQ-states for all different colored pairs. This has the effect, that 
equally colored quarks interact always repulsively, but different colors may attract or 
repell each other. An additional short-distance cutoff in the interaction is introduced: 
at r < 0.1 fm the potential is taken as a linear function of r, which results in a smooth 
quadratic decrease of the inter-particle force for r — > (soft-core model). This arbitrary 
cutoff can be justified in view of the finite spatial extension of the quark wave-function. 
The bookkeeping of all possible QQ pair- 



Table 2: Interaction table example for a 9-quark 
system. Dark and light gray fields are excluded due 
to the neglected self interaction and double counting, 
respectively. 



interactions in the system is realized using 
the so called "interaction table" of which 
an example is shown in table 121 The par- 
ticles are assigned a label and all particles 
have an associated color value. The ele- 
ments of the table show the (A* ■ Xj) values 
of the actual QQ pair realizations. The 
values of the D entries are randomly re- 
distributed periodically (see below). The 
total system is color-balanced (white), but 
neglecting the self-interaction results in the 
appearance of a net attraction. 

In the molecular dynamics simula- 
tion the quarks interact with the color- 
dependent effective pair-potential given by 
(JTJ. Periodic boundary conditions are ap- 
plied together with the Ewald summation 
technique for the proper treatment of the long range interaction. Forces acting on each of 
the particles (quarks) due to all other particles are calculated in every timestep and the 
equation of motion for each particle is integrated in time [|3]. 
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At the start of the simulation the particles are randomly distributed in the cubic sim- 
ulation box with kinetic energies determined by an initial temperature. An initialization 
phase is used before any measurement is taken on the system. In this phase the mo- 
menta of the particles are scaled back according to the given initial temperature in every 
timestep. The simulation time of this phase is chosen to be long enough to reach equilib- 
rium conditions for the unperturbed system over the subsequent measurement period. 

The effect of the gluon field on the quark component is taken into account in three 
ways: (i) it represents an infinite heat reservoir, energy can be coupled out of it without 
any changes in the qluon field; (ii) gluons can transfer color charge, the color state of two 
quarks may be exchanged (color rotation) with a characteristic time given by tq] (iii) qlu- 
ons perturb QQ-pairs, leading to a new realizations of the pair-interaction (redistribution 
of the interaction table) with a characteristic time td- 

The input parameters in our simulations are: density n = 10 quarks/fm 3 , initial tem- 
perature T = 200 MeV, effective quark mass m q = 300 MeV/c 2 , interaction coupling 
strength a = 1, time between reassignment of all D (different-color pair elements) in the 
interaction table td, time between color rotation events tq and the ratio of pairs color ro- 
tated in one event Nq- The plasma frequency is derived as up = 47ran/m q . The behavior 
of the system is now governed by the two, independently chosen "refreshment" times td 

Figure ^ shows the time- 
evolution of the quark kinetic 
temperature for different td val- 
ues. It can be seen, that in 
the limiting cases td = oo and 
td — ► the temperature fluctu- 
ates closely to its initial value. In- 
termediate values of td result in 
increasing temperatures. 

The r^rdependence of the tem- 
perature is plotted in figure E^a) 
for three snapshots in time. A 
remarkable resonant-like behavior 
of the heating of the system can 
be observed, where the maximum 
heating-rate establishes at td val- 
ues around the inverse plasma 
frequency 1/up ~ 0.19 fm/c. 
Figure 12(b) shows the dependence of the resonant-like peak on an additional color ro- 
tation rate, where tc = Td- We have found that the color exchange increases the energy 
transfer from the background field to the quark component by 10 to 20%. 

Figure E(a) displays the pair-correlation function (PCF) for the td = 0.0005 fm/c 
case neglecting color-rotation. The total PCF is decomposed into contributions of equal- 
colored-, and different-colored pairs. The attractive force acting between some of the 
different-colored pairs result in the appearance of the correlation peak at r pa 0. This 
indicates the onset of clusterization [ |H] of particles which gets more pronounced at 
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Figure 1 : Time-evolution of the kinetic temperature for sys- 
tems with different to values at tq = 0. (Numbers labeling the 
curves are the corresponding td values in fm/c units.) 
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Figure 2. Temperature vs. td at (a) times t = 1, 2 and 4 fm/c without color rotation; and (b) at time 
< = 4 fm/c with tc = td for color-rotation ratios Nc = 0, 50 and 100%. 



higher td values. Focusing only on the (always) repulsive equal-colored pairs, the PCF 
can be compared with classical OCP results [see figure E^b)]. This comparison shows 
that the repulsive equal-colored component of the QQ plasma is structurally similar to 
classical OCP in the so called "gas" phase (at T values around unity), which indicates 
already strong non-ideality of the QQ-plasma investigated here. 




Figure 3. Pair-correlation function for the td = 0.0005 fm/c case (a) decomposed into equal-colored 
and different-colored contributions; and (b) compared to classical OCP plasmas. Here aws = 0.3628 fm 
is the Wigner-Seitz radius, and T is the plasma (non-ideality) parameter. 



REFERENCES 

1. M. Gyulassy and L.McLerran, Nucl. Phys. A 750, 30 (2005). 

2. P. Levai, U. Heinz, Phys. Rev. C 57, 1879 (1998). 

3. D. Frenkel and B. Smit, Understanding Molecular Dynamics Simulations (Academic 
Press, New York, 2001). 

4. P. Hartmann, G. J. Kalman, Z. Donko and K. Kutasi, Phys. Rev. E 72, 026409 (2005). 

5. P. Hartmann, Z. Donko, G. J. Kalman and P. Levai, to be published. 

6. M. Thoma, these proceedings. 

7. E. V. Shuryak and I. Zahed, Phys. Rev. C 70, 021901 (R) (2004). 

8. M. E. Fisher and Y. Levin, Phys. Rev. Lett. 71, 3826 (1993). 



